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Abstract 

This paper aims at investigating the achievable performance and the issues that arise in ultra-dense 
networks (UDNs), when the signal propagation includes both the Line-of-Sight (LOS) and Non-Line- 
Of-Sight (NLOS) components. Backed by an analytical stochastic geometry-based model, we study 
the coverage, the Area Spectral Efficiency (ASE) and the energy efficiency of UDNs with LOS/NLOS 
propagation. We show that when LOS/NLOS propagation components are accounted for, the network 
suffers from low coverage and the ASE gain is lower than linear at high base station densities. However, 
this performance drop can partially be attenuated by means of frequency reuse, which is shown to 
improve the ASE vs coverage trade-off of cell densification, provided that we have a degree of freedom 
on the density of cells. In addition, from an energy efficiency standpoint, cell densification is shown 
to be inefficient when both LOS and NLOS components are taken into account. Overall, based on the 
findings of our work that assumes a more advanced system model compared to the current state-of-the- 
art, we claim that highly crowded environments of users represent the worst case scenario for ultra-dense 
networks. Namely, these are likely to face serious issues in terms of limited coverage. 

Index Terms 

Ultra-dense, LOS/NLOS, Area Spectral Efficiency, partially loaded, energy efficiency, coverage. 

I. Introduction 

There is a common and widely shared vision that next generation wireless networks will 
witness the proliferation of small-cells. As a matter of fact, researchers foresee network densifi¬ 
cation as one of key enablers of the 5' th generation (5G) wireless networks [[Q, [13. Although it 
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refers to a concept rather than being a precise definition, the term ultra-dense networks is used 
to describe networks characterized by a massive and dense deployment of small-cells, in which 
the amount of base stations may grow up to a point where it will exceed the amount of user 
devices 0. 

As wireless networks evolve, performance requirements for the new technology are becoming 
more and more stringent. In fact, the 5G requirements are set to a data rate increase up to a 
1000-fold with respect to current 4G systems [0, as well as for high energy efficiency 0 in 
order to limit the energy expenditure of network operators. Supported by recent results [j2l, 0, 
cell densification as been put forward as the main enabler to achieve these target data rates. For 
example, the authors in 0 have shown that the throughput gain is expected to grow linearly 
with the density of base stations per area; this is a result of the simplified system model used 
during the investigation. Namely, the assumption of a single slope path-loss model and that all 
base stations in the network are active. 

Nonetheless, further work on cell densification has shown that, under less ideal assumptions, 
network performance may be different than what predicted in 0. In particular, when different 
path-loss models than single slope are used, the actual performance of cell densification is less 
optimistic than what estimated with single slope path-loss 0-0. In addition, if the base station 
density increases beyond the user density, like in a typical ultra-dense scenario, the network has 
been shown to experience a coverage improvement at the expense of a limited throughput gain 
0, una. This implies that a larger number of BSs will need to be deployed to meet a given 
data rate target, translating on higher network infrastructures costs. 

A. Related Work 

In recent years, stochastic geometry has been gradually accepted as a mathematical tool for 
performance assessment of wireless networks. In fact, one of the most important contribution 
to the study of cell densification can be found in 0, where the authors proposed a stochastic 
geometry-based framework to model single-tier cellular wireless networks; by assuming a single 
slope path loss model, the authors have observed the independence of the Signal-to-Interference- 
plus-Noise-Ratio (SINR) and Spectral Efficiency (SE) from the BS deployment density, with the 
main consequence being the linear dependence of the ASE on the cell density. 

Nevertheless, when the assumption of single-slope path-loss is dropped, it emerges that SINR, 
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ASE and coverage exhibit a different behaviour than what was found in [|3]|. The authors 
in (6), lITTIl considered propagation models for millimeter waves. In fl6]|, the authors extended 
the stochastic geometry framework proposed in 0 to a multi-slope path loss model. The 
authors in lUTl developed a stochastic geometry framework for path-loss including Line-of-Sight 
(LOS) and Non-Line-of-Sight (NLOS) propagation for millimeter-waves. The effect of NLOS 
propagation on the outage probability has been studied in dTJ, where the authors propose a 
function that gives the probability to have LOS at a given point depending on the distance from 
the source, on the average size of the buildings and on the density of the buildings per area. A 
stochastic geometry-based framework to study the performance of the network with a combined 
LOS/NLOS propagation for micro-waves can be found in some previous work of ours flU] and 
in 0|. The common picture that emerges from the work is that of a non-linear behavior 

of the ASE with the cell density; overall, as a consequence of a different propagation model 
than the single slope path-loss, the spectral efficiency and coverage of the network do actually 
depend on the base station density. However, all these studies are based on the assumption of 
fully loaded networks, i.e., all the base stations are active and have at least one user to serve; thus, 
the applicability of the work above is limited to networks in which the density of base stations is 
lower than that of the users. In our paper, we broaden the study of network densification to the 
the case where there is no such a constraint in terms of BS density, i.e., we also tackle partially 
loaded networks, in which some base stations might be inactive. 

Work on stochastic geometry for partially loaded networks has been advanced in 0, IfTOll . 
EM . m. The authors in lfl2ll studied the coverage in single tier networks, while multi-tier 
networks are addressed in EM- An analysis of the area spectral efficiency of partially loaded 
networks has been carried out in [j3|j, while in IfTOll the authors have extended the stochastic 
geometry-based model further to include multi-antenna transmission, and have also assessed the 
energy efficiency. Overall, the authors in iflOl llT2ll . Ifl3ll have shown that the network coverage 
improves as the base station density increases beyond the user density; this is paid in terms of 
a lower throughput gain, which grows as a logarithmic function of the cell density. Nonetheless, 
the authors in IfTOll . Ifl2l . |[T3ll modeled the propagation according to a single slope path-loss 
model and did not investigate the effect of LOS/NLOS propagation in partially loaded networks. 

To the best of our knowledge, currently there is no work that addresses ultra-dense scenarios 
with path-loss models different than the single slope for both fully and partially loaded networks. 
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As a result of the combined effect of the path-loss model and of the partial load in ultra-dense 
networks, until now the behavior of the Area Spectral Efficiency (ASE), coverage, and energy 
efficiency as the networks become denser, was unknown. 

B. Our Contribution 

This paper seeks to investigate the cell densification process in ultra-dense networks and 
evaluate the effect of LOS/NLOS propagation on performance metrics such as coverage, spec¬ 
tral efficiency, area spectral efficiency, and energy efficiency. Specifically, we use ultra-dense 
networks as a term to refer to those networks characterized by a very high density of base 
stations and that include both the cases of fully loaded networks (i.e., networks in which all 
the base stations are active) and partially loaded networks (i.e., networks in which some base 
stations might be inactive and not transmit to any user). Overall, the major contributions of our 
work can be summarized in the following points: 

1) Stochastic geometry-based model for ultra-dense networks with LOS/ NLOS propa¬ 
gation: We propose a model based on stochastic geometry that allows us to study the Signal-to- 
Interference-plus-Noise-Ratio (SINR) distribution, the spectral efficiency and the area spectral 
efficiency of ultra-dense networks where the propagation has LOS and NLOS components. We 
build on previous work |0 and we adapt the model proposed by Andrews et al. to the case of 
LOS/NLOS propagation. In addition, our framework takes also into account the partial load of 
ultra-dense networks, in which a fraction of base stations may be inactive. 

2) Study of cell densification, partial load and frequency reuse in networks where signal 
follows LOS/NLOS propagation: Lirst, we investigate the effect of network densification on 
performance metrics such as SINR, spectral efficiency and area spectral efficiency in networks 
where the path-loss follows the LOS/NLOS propagation. In particular, the ASE gain becomes 
lower than linear at high cell densities, meaning that a larger number of BSs would be necessary 
to achieve a given target with respect to the case of single slope path loss. Moreover, the network 
coverage drops drastically as the BS density increases. Then, we show that the performance drop 
due to LOS/NLOS propagation gets mitigated by the usage of frequency reuse or if the base 
station density exceeds the user density, as it is likely to occur in ultra-dense networks. To the 
best of our knowledge, the combined effect of LOS/NLOS propagation and partial load/frequency 
reuse has not been addressed before. 
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3) Investigation on the minimum transmit power per BS and energy efficiency for 
networks where signal propagates according to LOS/NLOS path-loss: As the cell density 
increases, the transmit (TX) power per base station can be lowered. We evaluate the minimum 
TX power per BS such that the network is guaranteed to be in the interference-limited regime, 
in which case the performance is not limited by the TX power. Second, we make use of the 
TX power to determine the energy efficiency of the network when the propagation has LOS 
and NLOS components. We show that the energy efficiency with LOS/NLOS propagation drops 
considerably with respect to the case of single slope path-loss, making cell-densification costly for 
the network from an energetic stand-point. We further extend the study of the energy efficiency 
to frequency reuse and partial load. 

C. Paper Structure 

The remainder of this paper is organized as follows. In Section [TT] we describe the system 
model. We show our formulation for computing the SINR, SE and ASE in Section HIT] and we 
address the energy efficiency in Section [TV] In Section [V] we present and discuss the results 
while the conclusions are drawn in Section [VH 


II. System model 


In this paper we consider a network of small-cell base stations deployed according to a 
homogeneous and isotropic Spatial Poisson Point Process (SPPP), denoted as <f> C M 2 , with 
intensity A. Further, we assume that each Base Station (BS) transmits with an isotropic antenna 
and with the same power, Ptx, of which the value is not specified, in order to keep our model 
general and valid for different base station classes (e.g., micro-BSs, pico-BSs, femto-BSs); we 
focus our analysis on the downlink. 

A. Channel model 

In our analysis, we considered the following path loss model: 



K\d /3l with probability p L (d), 
Kxhd~^ with probability 1 — Ph{d), 


(1) 
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where /? L and /3 NL are the path-loss exponents for LOS and NLOS propagation, respectively; 
K\ i and A" NL are the signal attenuations at distance d — 1 for LOS and NLOS propagation^ 
respectively; p L (cf) is the probability of having LOS as a function of the distance d. The model 
given in (Q]) is used by the 3GPP to model the LOS/NLOS propagation, for example, in scenarios 
with Heterogeneous Networks lfT4l Table A.2.1.1.2-3]. The incorporation of the NLOS component 
in the path loss model accounts for possible obstructions of the signal due to large scale objects 
(e.g. buildings ), which will result in a higher attenuation of the NLOS propagation compared 
to the LOS path. We further assume that the propagation is affected by Rayleigh fading, which 
is exponentially distributed ~ exp(/i). 

Regarding the shadow fading, it has been shown that in networks with a deterministic, either 
regular or irregular, base station distribution affected by log-normal shadow fading, the statistic 
of the propagation coefficients converges to that of a network with SPPP distribution as the 
shadowing variance increases lfl5l . In other words, this SPPP intrinsically models the effect of 
shadow fading. 

B. LOS probability function 

To ensure that our formulation and the outcomes of our study are general and not limited to a 
specific LOS probability pattern, we consider two different LOS probability functions. The first 
one, which is proposed by the 3GPP Ifl4l Table A.2.1.1.2-3] to assess the network performance 
in pico-cell scenarios, is given below: 

Pl,3g(^0 = 0.5 — min (0.5,5e“~^) + min (0.5,56”^), (2) 

with d 0 and d\ being two parameters that allow d2]) to match the measurement data. Unfortunately, 
this function is not practical for an analytical formulation. Therefore, we chose to approximate 
it with a more tractable one, namely: 

p L (d) = exp (— (d/L) 2 ) , (3) 

where L is a parameter that allows ([3]) to be tuned to match ©. The second function is also 
suggested by the 3GPP HU Table A.2.1.1.2-3] and is given below: 

Ph(d) = exp (-d/L). (4) 


'The parameters Kl and TTnl can either refer to the signal attenuations at distance d = lmord=l km; this depends on 
the actual values given for the parameters of the channel model. 
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From a physical stand point, the parameter L can be interpreted as the LOS likelihood of a given 
propagation environment as a function of the distance. 

C. User distribution, fully and partially loaded networks 

In our model, we always assume that: (i) the users are uniformly distributed; (ii) that each 
user’s position is independent of the other users’ position; and (iii) each user connects only to 
one base station, the one that provides the strongest signal. We denote by Au the density of 
users per area; whenever we consider a finite area A, N\j indicates the average number of users 
in the network. We also assume the users are served with full buffer, i.e., the base station has 
always data to transmit to the users and make full use of the available resources. 

Depending on the ratio between the density of users and the density of base stations, we dis¬ 
tinguish between two cases, namely, fully loaded and partially (or fractionally) loaded networks. 
By fully loaded networks we refer to the case where each BS has at least one user to serve. 
With reference to a real scenario, fully loaded networks model the case where there are many 
more users than base stations, so that each base station serves a non-empty set of users. 

However, when the density of users is comparable with or even less that of base stations, 
some base stations may not have any users to serve and would then be inactive, meaning that 
they do not transmit and do not generate interference. When this occurs, we say that the network 
is partially loaded. The network can be modelled as partially loaded to study those scenarios 
characterized by high density of base stations and, in particular, scenarios where the density of 
base stations exceeds the density of users, such as in ultra-dense networks. 

To define formally the concepts of fully and partially loaded networks, we first need to 
introduce another concept, which is the probability a base station being active. 

Definition 1 (Probability of a base station being active). The probability of a base station being 
active, denoted as p&, is the probability that a base station has at least one user to serx’e. This 
event implies that the base station is active and transmits to its users. 

Definition 2 (Fully loaded and partially loaded networks). The network is said to be fully loaded 
tf Pa = 1; the network is said to be partially or fractionally loaded if Pa < 1. 


III. SINR, SPECTRAL EFFICIENCY AND ASE 

In this section we propose an analytical model to compute the SINR Complementary Cumu¬ 
lative Distribution Function (CCDF), which allows us to asses key performance metrics such as 
coverage, spectral efficiency and ASE. 


A. Procedure to compute the SINR CCDF 

In order to compute the SINR tail distribution (i.e., the Complementary CDF), we extend the 
analytical framework first proposed in [Q so that to include the EOS and NFOS components. 
From the Slivnyak’s Theorem [fT6l Theorem8.10], we consider the typical user as the focus of 
our analysis, which for convenience is assumed to be located at the origin. The procedure is 
composed of two steps: (i) we compute the SINR CCDF for the typical user conditioned on 
the distance from the user to the serving base station, denoted as r; (ii) using the PDF of the 
distance from the closest BS f r (R), which corresponds to the serving BS, we can average the 
SINR CCDF over all possible values of distance r. 

Fet us denote the SINR by 7 ; formally, the CCDF of 7 is computed as: 

f+O o 

p [7 > y] = E r [P [7 > y\r\ ]= P['Y>y\r = R] f r (R)dR. (5) 

Jo 

The key elements of this procedure are the PDF of the distance to the nearest base station f r (R) 
and the tail probability of the SINR conditioned on r, P [7 > y\r = R\. The methodology to 
compute each of these elements while modelling the FOS and NFOS path loss components will 
be exposed next. 


B. SPPPs of base stations in LOS and in NLOS with the user 

The set of the base stations locations originates an SPPP, which we denote by <3> = {x n }o As 
a result of the propagation model we have adopted in our analysis (see Section III-AI) . the user 
can either be in FOS or NFOS with any base station x n of <I>. Now, we perform the following 
mapping: we first define the set of FOS points, namely <f> L , and the set of NFOS points, <f> NL . 
Then, each point x n of $ is mapped into if the base station at location x n is in FOS with 
the user, while it is mapped to <3> NL if the base station at location x n is in NFOS with the user. 
Since the probability that x n is in FOS with the user is p L (IMI)> it follows that each point x n 


2 Whenever there is no chance of confusion, we drop the subscript n and use x and instead of x n for convenience of notation. 
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of $ is mapped with probability p L (||x||) into <f> L and probability Pnl(IMI) = 1 — Pl(IMI) into 
$ NL . Given that this mapping is performed independently for each point in <I>, then from the 
"Thinning Theorem" Ifl6l Theorem 2.36] it follows that the processes <f> L and <f> NL are SPPPs 
with density A L (x) = Ap L (IMI) and A NL (x) = A (1 — p L (||x||)), respectively. Note that, because 
of the dependence of A L (x) and A NL (a;) on x, <f>L and $ NL are inhomogeneous SPPPs. To make 
the formulation more tractable, we consider $ L and <3> NL to be independent processes; because 
the union of two SPPPs processes is an SPPP of which the density is the sum of the densities of 
the individual SPPPs iflTl Preposition 1.3.3], the union of <3> L and <f> NL is an SPPP with density 
Al(x) + Anl(^) = A, i.e., it is an SPPP with the same density as that of the original process <3>. 
Hence, the assumption of independence between and $nl does not alter the nature of the 
process <f>. 

C. Mapping the NLOS SPPP into an equivalent LOS SPPP 

Given that we have two inhomogeneous SPPP processes, it is not trivial to obtain the distri¬ 
bution of the minimum distance of the user to the serving base station, which will be necessary 
later on to compute the SINR CDF. In fact, assuming the user be in LOS with the serving base 
station at a distance d\, there might be an interfering BS at a distance d 2 < d\ which is in NLOS 
with the user. This is possible because the NLOS propagation is affected by a higher attenuation 
than the LOS propagation. 

Hence, to make our problem more tractable, we map the set of points of the process $nl, 
which corresponds to the NLOS base stations, into an equivalent LOS process $ eq ; each point 
x 6 $nl located at distance c/nl from the user is mapped to a point x eq located at distance d eq 
from the user, so that the BS located at x eq provides the same signal power to the user with 
path-loss K L d~P L as if the base station were located at x with path-loss Ji N l^nl NL - 

Definition 3 (Mapping function f eq ). We define the mapping function f eq : <f> NL — y <f> eq as: 



x 


( 6 ) 



( 7 ) 
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Definition 4 (Inverse mapping function g eq ). The inverse function g eq = f cq : <f> eq —> $ NL is 
defined as: 

9eq(x) = M^^q 1 (IM|) , (8) 

11 CC 11 

<(d) = l^j d^/^, = Keqd 8~ l} (9) 

/ \ 1//3nl 

where K eq = while (3 eq = /3 l //3 N l- 

It is important to notice that, from the "Mapping Theorem" lfT6l Theorem 2.34], $ eq is still 
an SPPP. 

1). PDF of the distance from the user to the serving BS 

Using the mapping we introduced in Section UlI-Cl we can compute the PDF f r (R) of the 
minimum distance r between the user and the serving BS. We first compute the probability 
P [r > 5]; the PDF can be ultimately obtained from the derivative of P [r > R] as f r (T>) = 
^(1 — P [r > R]). Let /i(0, V) be the ball of radius l centred at the origin (0,0). Moreover, 
we use the notation $(^4) to refer to number of points x G $ contained in A lfl6ll . Using the 
mapping we introduced in Section UlI-Cl the probability P [r > 7i] can be found as: 

P [r > R] = P [<P L (B( 0, R))= OH <F eq (.B (0, R)) = 0] 

( = } P [<P L (5(0, R)) = on <P NL (B (0, <(5))) = 0] 

( = } P [<P L (5(0, R)) = 0] • P [<P NL (5 (0, <(5))) = 0] , (10) 

where equality (a) comes from the mapping defined in ([8]) and in ([9]), while equality ( b ) comes 
from the independence of the processes 'Ll an d $nl- By applying the probability function of 
inhomogeneous SPPP |jT6l Definition 2.10], we obtain the following, 

P [r > 5] = exp ( — I Al(x)Tx J exp ( — ANiXaOdx ). (11) 

V Jb{o,r) J \ JB(o,d^( R )) / 

From (fTTh . we can obtain /,.(/?), first, by integrating and, second, by computing its first derivative 

in R. The formulation in (fTTh is general and thus can be applied to several LOS probability 

functions Ph(d). Below, we provide the expression of the PDF of the distance from the UE to 

the serving BS for the LOS functions ([3]) and ©, respectively. 
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Result 1. If the LOS probability function is as in © and if we denote d eq (R) by R cq , the PDF 
of the distance to the serving BS is: 


fr(R) 


tA L 2 e 


R 

'77 


—7rAL 2 e L 2 


e ~-K\Rl q 


( 12 ) 


-K^R 2 ^ 


—2nXRe~^ 7rAA"g q 2 / 3 eq i? 2/3eq_1 e-^-- nXK^^R 2 ^- 1 ) . 


q ” e q 


Result 2. If the LOS probability function is as in © and if we denote d cq (R) by R cq , the PDF 
of the distance to the serving BS is: 


fr(R) = - e 


^2irXL 2 e Jl^XLRe —'kXR* ~—2ttXL 2 e 

• c *6 eq • e 


. g—2-7TA LR eq e 


(13) 


( - 2ttA Le~% - 2vrA(L - R)e~% - nXK^^R 2 ^- 1 

„ KeqRPen K eq RP«<l \ 

+ 2 TrXLK cq P eq R^e -- + 2 nXLK eq f 3 eq R^- 1 (K eq R^ - L)e - —J. 

We refer to the Appendix for the details of the f r (R) we have given in (fl2l) and in (fl3l) . 


E. Spatial process of the interfering of the active base stations 

In this framework, we include the cases of partially loaded networks and of frequency reuse 
and we treat them separately. First, we identify the set of BSs that are active, i.e., those having 
one or more users to serve. As we focus the analysis on the typical user, we can also identify 
the set of BSs that act as interferes for that user; an active BS (excluding the one serving the 
user) is seen as an interferer if that BS transmits over the same band used to serve that user. 
In the following, we denote by <1>a the set of active BSs, while we denote by <f>i the set of the 
interfering BSs. 

Let us consider first the case of frequency reuse, in which all the base stations are active, but 
each of these only uses a portion of the spectrum, in order to reduce the interference. Since 
all the base stations are supposed to be active, the process $a is the same as $. However, we 
assume each base station selects a channel in a random manner using, for instance, frequency 
ALOHA spectrum access lfT8ll . With a frequency reuse factor of N, each base station uses 1 
out of N channels, which is chosen independently from the other base stations. Hence, each BS 
interferes with a given user with probability 1/iV; this is equivalent to carrying out a thinning 
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of the original process <3> with probability 1/N; from the Thinning Theorem, we obtain that <1>\ 
is a homogeneous process with density Ai = X/N. 


In regards to the partially loaded networks, we recall from Section III-CI that a fraction of the 


base stations might be inactive and, as such, would not generate interference. Assuming all the 
BSs transmit over the same band, in partially loaded networks the active base stations are the 
only BSs that generate interference to the users—with the exception of the serving BS. Thus, we 
can write $i = <3 >a \ %o, where xq is the serving base station; moreover, from the Palm Theorem 
m, and $a are characterized by the same density. To obtain the process of active BSs <1> A 
from the original the process <1\ we first assume that each UE deployed in the network connects 
to the closest BS; 3 finally, only the BSs which are assigned one or more users are active. Yet, 
this is equivalent to performing a thinning of the original process to obtain <h A . However, the 
fact that a base station is picked to be part of <3 >a (i.e. there exists a user for which this BS is 
the closest one) depends on the positions of the neighbouring base stations, which implies that 
the base stations are not chosen independently of one another Ifl2l . As the independence is one 
of the necessary conditions in order to have an SPPP, it follows that $ A is not an SPPP 
Although <3 >a cannot be formally regarded as an SPPP, it has been shown that the actual SPPP 
obtained through the thinning of $ well approximates <h A IflOl . |[T2l . Specifically, the authors in 
lfl2l have shown that, (i) the probability p A of a base station to be active (i.e., to have users to 
serve) can be well approximated once the density of users Au and density of base stations A are 
known, and, (ii) the process <h A of active base stations can be well approximated by thinning 
the original process <1 >a with probability p A , which is given below l!T2ll : 



(14) 


From the Thinning Theorem, it follows that the resulting process obtained through thinning as 
described above is an SPPP with density Aa = PaA; moreover, Ai = Aa- In light of these 
findings, we approximate $a with an SPPP of density Aa = PaX. 

Fig. H] shows how the probability p A and the intensity of the interfering BSs Ai vary as 
functions of the ratio A/Au- 

3 As we recall from Section IIII-CI with LOS/NLOS propagation the serving BS might not be the closest one to the user. 


Nonetheless, this does not alter the validity of the explanation we are giving in this section. 
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Fig. 1. Probability of a BS being active and density of interfering BS vs BS density for partially loaded networks. The 
probability pa drops as the ratio A/Au is close to or greater than 1, i.e., as A approaches Au- As a result of this, the density 
of active BSs as well as the density of interfering BSs converge to Au as A approaches or overcomes Au- 


F. SINR inverse cumulative distribution function 

The probability P [7 > y\r = R] can be computed as in [[51 Theorem 1]; we skip the details 
and provide the general formulation: 

~gK L R-fo 


p [7 > y\ r — R] — p 


> y 


_ -pyK^R^a 2 


= e 


C^yKpR^), 


(15) 


<T 2 + Ir 

where g is the Rayleigh fading, which we assume to be an exponential random variable ~ expi/i); 
a 2 is the variance of the additive white Gaussian noise normalized with the respect to the transmit 
power, I R is the interference conditioned on R, i.e., 


Ir= Y 9i K A\ x i\\ ^ + Y 9j K L\\xj\\ /3l (16) 

{i: 31ie < f>Ln<f>A, |pi||>H} {j: /eq(31j)G3>NLnT>A, ll/eq(aij)||>K} 

= Y 9iKi]\xi\\-^ + Y ^iTNL||^||-^ NL (17) 

{ i : XiS'f’LnT-A, ||a:i||>/?} {j : i ;; G<f>NLn$A, \\xj I^^Jq 1 (R)} 

where gi and g 3 are independent and identically distributed ~ expi/i) fading coefficients. Please, 
note that the overall interference accounts only for the active base stations. Compared to the 
formulation of £i r (s) proposed in 0, in our case we have to deal with two non-homogeneous 
SPPP, namely <f> L and $ NL instead of a single homogeneous SPPP. The Laplace transform C Ir (s) 
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can be written as follows: 


£-Ir( S ) E'^Ln$A,'I > NLn$A,3i,3j 


exp | -s ^ giK L \\xi\\ ^ 

{i\ lig^Ln^A, |®t|>l®o|} 


exp ( — s ^ gj K NL ll^ll /5nl 

{j: ij€$NLn$A, ||ajj||>rf^-(-R)} 

Given that <f> L and <f> NL are two independent SPPP, we can separate the expectation to obtain: 


^g( s ) — E$ L n$ A ,9i 


exp l — s ^2 9i K h\\xi\\ ^ 

{*: lie^LO'I'A, ||a:i||>fl} 

I —/?nl 


(18) 


E $NLn^A,Si exp|^-s ^2 9jKNL\\xj\ 

{j:xj£® nlH$a, \\xj (-R)} 

By applying the Probability Generating Functional (PGFL) for SPPP (which holds also in case 
of inhomogeneous SPPP Ifl6l ) to (fl8l) . we obtain the following result: 

Result 3. The Laplace transform Ci R (s) for LOS/NLOS propagation with model given in (]T|) is: 

+oo 


£i r (s) = exp - 2n\\ 


R 


sK]V 


-Pc 


sKi/o Pl -f n 


Ph(v)vdv 


+00 


exp — 2n\ 


def(R) 


sKxlV ^ 

sKf^i,V~P Nh + /i 


p NL (v)vdv 


(19) 


The Laplace transform in (fl9l) along with tlTTb and (IA.4I) can be plugged in © to obtain the 
SINR CCDF through numerical integration. 


G. Average Spectral Efficiency and Area Spectral Efficiency 

Similarly to © Section IV] we compute the average spectral efficiency and the ASE of the 
network. First, we define the ASE as: 

Aa • A • BWu • E[C] A a • E[C] 

= — jrm: — = —iv—■ (20) 

where BWa is the available bandwidth, BWu is the used bandwidth, E[C] is the average spectral 
efficiency, A is the area, Nbs,a is the number of active base stations within the area A, and N 
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is frequency reuse factor. The average rate E[C] can be computed as follows 0: 


r»+oo 


E[C] = E [log 2 (l + 7 )] = 

/*+oo /*+oo 


P [log 2 (l + 7 ) > u] d u 


( 21 ) 


'0 JO 


r»+oo /*+oo 


'0 JO 


P [log 2 (l + 7 ) > u\r = R] f r (R)dRdu. 
e - M ( 2 --i c Ir {h(2 u - l)K~ l R^)f r (R)dRdu 


( 22 ) 


where jCi r {s) is given in (fl9l) . Similarly to the SINR CCDF, (1221) can be evaluated numerically. 


IV. Energy efficiency with LOS/NLOS propagation 
A. Computing the transmit power per base station 

We evaluate the TX power in order to compute the overall power consumption of the wireless 
network. Ideally, to ensure the network performance not be limited by the transmit power, Ptx 
should be set in order to guarantee the interference-limited regime, i.e., the transmit power should 
be high enough so that the thermal noise power at the user receiver can be neglected with respect 
to the interference power at the receiver. In fact, when the network is in the interference-limited 
regime, the transmit power is high enough that any further increase of it would be pointless in 
terms of enhancing the SINR, since the receive power increment would be balanced by the exact 
same interference increment. 

Practically, we refer to the outage probability 9 = P [7 < 7 th] as a constraint to set the power 
necessary to reach the interference limited regime. When the TX power is low, small increments 
of Ptx yields large improvements of the outage 9\ however, as Ptx increases, the corresponding 
outage gain reduces, until 6 eventually converges to its optimal value 9*, which would be reached 
in absence of thermal noise. It is reasonable to assume the network be in the interference-limited 
regime when the following condition is met: 

1 9* -9\< A e , (23) 

where Aq is a tolerance measure setting the constraint in terms of the maximum deviation of 9 
from the optimal value 9*. Eq. (1231) provides us with a metric to compute the transmit power, but 
does not give us any indication on how to find P T x as a function of the density A. Unfortunately, 
we cannot derive a closed-form expression for the transmit power that satisfies (1231) . as we do 
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not have any closed-form solution for the outage probability 6 = P [7 < 7 th ], We then take a 
different approach to calculate the minimum transmit power. 

In Alg. [I] we proposed a simple iterative algorithm that finds the minimum transmit power 
satisfying (l23l) by using numerical integration of ©. This algorithm computes the outage 
probability corresponding to a given P TX ; starting from a low value of power, it gradually 
increases Ptx by a given power step Ap, until (l23l) is satisfied. To speed up this procedure, the 
step granularity is adjusted from a coarse step pi up to the finest step p X[) , which represents the 
precision of the power value returned by Alg. [0 


Algorithm 1 Steps to compute the transmit power. 

INPUTS: 

1) Vector of the power steps in dBm p = [pi, ■ • • pat p ], N p is the length of vector p; 

2) Outage SINR threshold 7 th ; 

3) Outage tolerance A g\ 

Initialize variables: 

1) Pcurr — Pn 0 , where P\r 0 is the AWGN power in dBm over the bandwidth BWu 

2 ) Pfin Pcurr 

Find optimal outage 9* = P [7 < 7 th] by integrating 0 with parameter a 2 = 0 
for k = 1, ■ • ■ , N p do 

Find 9(P curr ) = P [7 < 7 th] by integrating (0 with parameter a 2 = 10 
Set granularity of the power step A p = p/, 

while \d* - 9{P curr ) | > Ag do 

Increase the current power with step Ap, i.e, P curr = Pcurr + Ap 
Find 9(P curr ) = P [7 < 7 t J by integrating 0 with parameter er 2 = 10 
Update the final value of power, i.e., Pa n = P curr 
Remove the last power increment before increasing the granularity, i.e., P curr = P CU rr — Ap 
OUTPUT: Pfi n is the power in dBm s.t. (l23l) is satisfied. 


B. Energy efficiency 

In this subsection, we characterize the energy efficiency of the network as a function of the 
base station density to identify the trade-off between the area spectral efficiency and the power 
consumed by network. We define the energy efficiency as the ratio between the overall throughput 
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delivered by the network and the total power consumed by the wireless network, i.e., we define 
the energy efficiency as follows: 


Pee (A) — 


T( A) 


(24) 


Ptot(A) ’ 

where T(A) is the network throughput, which can be written as T(A) = A ■ BW • r/ A {X), with 
BW denoting the bandwidth and //a(A) denoting the area spectral efficiency; Ptot is the overall 
power consumption of the network. 

When we compute the power consumption of each BS, we need to take into account that 
a fraction of the base stations may be inactive and model the power consumption accordingly. 
For active base station, we model the power consumption Pbs,a of the base station assuming 
that Pbs,a is the sum of two components, i.e., Pbs,a = Po + Prf; the first, denoted by P 0 , 
takes into account the energy necessary for signal processing and to power up the base station 
circuitry. This power P 0 is modelled as a component being independent of the transmit power 
and of the base station load m. The second component, denoted by P RF , takes into account 
the power fed into the power amplifier which is then radiated for signal transmission. The power 
Prf is considered to be proportional to the power transmitted by the base station; we can thus 
write Prf = P rf Ptx> where P RF takes into account the losses of the power amplifier (i.e., we 
assume Prf to be the inverse of the power amplifier efficiency). 

In the case of inactive base stations, we assume that the BS switches to a stand-by state 
for energy saving purposes Il20ll . in which it does not transmit (i.e., P RF = 0) and reduces the 
circuitry power consumption. Therefore, the power required to maintain the stand-by state can 
be modelled as Pbs,s = pPo, where p is power saving factor that describes the relative power 
consumption of the circuitry with respect to the active case; note that 0 < p < 1- 

Finally, the overall power consumption of the network due to both active and inactive base 
station can be expressed as follows: 


Ptot = AX a Pbs,a + Al(A — Aa)Pbs,s 
= AX a P 0 + AAaPtxPrf + Al(A — X A )pP 0 . (25) 


The energy efficiency for the cases of fully loaded networks and partially loaded networks is 
addressed in the following sub-sections. 
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C. Energy efficiency for fully loaded networks 

In this section we study the energy efficiency r] EE (X) trend as a function of A; we focus on 
fully loaded networks, i.e., j>a = 1 and A a = A. Unfortunately, the analysis of the derivative of 
t]ee is not straightforward, as we have a closed-form solution neither for the throughput T(A) 
nor for the transmit power Prx(A). One feasible way to get around this burden is to approximate 
T( A) and Prx(A) with functions in the form: 


f(z) = az b . 


(26) 


The model in (1261) has two advantages: (i) it can be easily derivated and, thus, is apt to investigate 
the existence of optima; (ii) it is well suited to fit the non-linear behaviour of ASE and TX power. 
In fact, we have shown in our previous work [f2Tll that both T(A) and Prx(A) can be approximated 
with a piece-wise function in the form (l26l) . where the parameters a and b can be obtained, for 
instance, by linear regression in the logarithmic domain for a given range of values of A. 

Backed by the conclusions from previous work [|2T|. we approximate the throughput as T(A) = 
AT 0 X a and the transmit power as Prx(A) = PrA 5 , within a given interval of A. Under these 
assumptions, the energy efficiency becomes: 

T 0 X a io' 


%e(A) — 


TnA"" 1 


(27) 


(28) 


APo + AArf-PtA * 5 Po + Arj'PtA ' 5 

The derivative of ??ee(A) is given below: 

cW(A) = T 0 P 0 (a - l)X a ~ 2 + T 0 K RF P T (a - 5 - l)X a+s ~ 2 
dA (P 0 + ArfPtA - 5 ) 2 

Let us note that T 0 , P 0 , A" RF and P T are positive; moreover it is reasonable to assume that a > 0 
(i.e., the area spectral efficiency is an increasing function of the density) and that 5 < 0 , i.e., 
the transmit power per BS is a decreasing function of the density. In the following paragraphs, 
we study the behaviour of the energy efficiency as function of the density A by analyzing the 
derivative r/ EE (A). We distinguish the following three cases. 

1) The energy efficiency is a monotonically increasing function: This occurs if the ASE growth 
is linear or superlinear, i.e., if a > 1. From this, if follows that also a > 1 > 1 + 5 holds true; in 
this case, r/^fA) is strictly positive, meaning that the energy efficiency increases as the density 


increases. 
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2) The energy efficiency is a monotonically decreasing function: This occurs if the ASE 
growth is sublinear, i.e., if a < 1, and, in addition, a < 1 + 5. Then, rj' EE ( A) is strictly negative 
and so the energy efficiency is a monotonically decreasing function of the density A. 

3) The energy efficiency exhibits an optimum point: If ASE gain is sublinear (i.e. a < 1) 
but grows with a slope a sufficiently high, (i.e., a > 1 + 5), then we obtain that the derivative 
^ee(A) nulls for 



(29) 


is positive for A < A 0 and is negative for A > A 0 . Therefore, A 0 is a global maximum of the 
energy efficiency. 


As a whole, the behavior of the spectral efficiency is due to how the growths of the ASE of 


the TX power relate among each other as the base station density increases. If the ASE grows 
rapidly enough to counterbalance the total power increase of the network given by the addition 
of base stations, then the energy efficiency increases with the BS density; this means that adding 
base station is profitable in terms of energy efficiency. Else, adding BSs turns not to be profitable 
from energy efficiency point of view. 

D. Energy efficiency for partially loaded networks 

For partially loaded networks, we only analyze the case where A > Au, as the opposite case 
of A < Au leads back to fully loaded networks. By using L’HopitaTs rule, one can show that 
(fl4l) can be approximated by p\ = AjjA - 1 , for A is sufficiently greater than Au- By applying 
this approximation to (1251) . we obtain: 


Ptot — Au-Pq(1 ~ P) + -V-Po + ^u-Prf-Pt-^ 5 - 


(30) 


It is known from Ifl9l that, as the BS density increases, the main contribution to the total power 
consumption is due to the circuitry power P 0 , while the transmit power becomes negligible 
for the overall power balance. Therefore, to make the problem more tractable, we can further 
approximate the total power in (1301) as / jot — AuPo(l — p) + ApP 0 . From (1241) . by using the 
approximation T(A) = AT 0 A" for the throughput and P TO t — AuPo(l — p) + ApP 0 for the power, 
we obtain the following expression for the energy efficiency: 


77EE(A) ~ AuPo(l-p) + ApP 0 ' 


(31) 
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To analyze the behaviour of the energy efficiency as a function of A, we follow the same approach 
as in Section HV-CI and we compute the derivative of p E e(A), which is given below: 

dr/EE^) _ TpX*- 1 (A p(a — 1) + aAu(l — p )) 

dA (AuP 0 (l - P) + ApP 0 ) 2 

As the ASE is known to be sub-linear for partially loaded networks 0, iflOll . we assume 0 < 
a < 1; moreover, the power saving factor p satisfies 0 < p < I- Therefore, the derivative rf EE 


nulls for: 


A* = 


aAu(l — p) 


(33) 


p(l-a) 

it is positive for A < A* while it is negative for A > A*. Hence, A* is a local maximum of the 
energy efficiency for partially loaded networks and the energy efficiency decreases for densities 
A > A*. Note that, this result holds for A sufficiently greater than Ay. 


V. Results 

In this section we present and discuss the results we obtained by integrating numerically the 
expressions of outage probability, of the Spectral Efficiency (SE), and of the ASE. In Section 
IV-A1IV-BI and IV-Cl we assume the network to be interference-limited (i.e., we set the thermal 
noise power to 0), while the noise is taken into account in Section IV-DI and IV-E1 
The parameters we used to obtain our results are specified in Table [D 

TABLE i 

Parameters for results section 


Parameter 

Value 

Path-loss - Single slope 

PL SL (d km ) = 140.7 + 36.71og(d km ), P = 3.67, JTsl = 10 1407 (H) 

Path-loss - Combined LOS/NLOS 

See Q: with d in km, AT = 10 1038 , /?l = 2.09, ATm = 10 14 54 , Pnl = 3.75, 
do = 0.156km, d\ = 0.03km Il4l 

Parameter L 

82.5m, set so that (2) and (3) intersect at the point corresponding to probability 0.5 

Bandwidth BW 

10 MHz centered at 2 GHz 

Noise 

Additive White Gaussian Noise (AWGN) with -174 dBm/Hz Power Spectral Density 

Noise Figure 

9 dB 

Antenna at BS and UE 

Omni-directional with 0 dBi gain 

A'rf 

10 Q9] 

Po 

10W (H 




























Outage probability - P[y^ Y th ] 
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A. Spectral efficiency, outage probability and ASE 

In this subsection we assume the network to be fully loaded and with frequency reuse 1. We 
compared the results for two LOS probability functions, namely © and ©; we also compared 
the results for LOS/NLOS propagation with those obtained with a the single slope path-loss 
model. We first analyze the outage probability (defined as 9 = P [7 < 7 th]) results, which have 
been obtained by numerical integration of ( 0 ). 




(a) Outage vs base station density. (b) ASE vs base station density. 

Fig. 2. (a) Outage probability and (b) ASE vs base station density for different LOS probability function. 

We show the outage probability results in Fig. [2al where we can see the different effect of the 
LOS/NLOS propagation with respect to the single slope Path-Loss (PL). With single-slope PL, 
the outage is constant with the BS density. In contrast, with LOS/NLOS propagation, there is a 
minimum in the outage curves, which is achieved for density A = 50-100BSs/km 2 , depending 
on the LOS probability function. Within this range of densities, the user is likely to be in LOS 
with the serving BS and in NLOS with most of the interfering BS, meaning that the interference 
power is lower than the received power. 

At densities A greater than 300BSs/km 2 , the outage starts growing drastically and, depending 
on the LOS likelihood, can reach 32-43%. This is due to the fact that more and more interfering 
BSs are likely to enter the LOS region, causing an overall interference growth and thus a reduction 
of the SIR. At densities A smaller than lOOBSs/km 2 , the serving BS as well as the interfering 
BSs are likely to be in NLOS with the user. Because of this, both the receive power and the 
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overall interference increase at the same pacqj and, as a consequence, the SIR remains constant, 
and so does the outage. Let us note that, the LOS probability function affects the outage curves 
at intermediate values of the BS density (e.g. 10-300 BSs/km 2 ). At low densities, all the BSs 
are likely to be in NLOS with the user, while at high densities the serving BS and the strongest 
interferes are likely to be in LOS BSs are likely to be in LOS with the user. 

The results of the ASE are shown in Fig. [Zb] Compared to the single-slope PL, which shows 
a linear growth of the ASE with the density A, with the LOS/NLOS propagation we observe a 
different behaviour of the ASE. In particular, we observe a lower steepness of the ASE curve 
at high BS densities, which is due to the effect of the interfering BSs entering the LOS region 
and, thus, increasing the total interference power. 

To assess steepness of the ASE, we can use linear regression to interpolate the ASE curve 
with the model given in (1261) . In particular, we can approximate the ASE t/ a (A) with a piece-wise 
function of the kind 7 a (A) = r/ AU \ n , where r/ A 0 and a are given for given intervals of A. We 
specifically focus on a, which gives the steepness of the ASE curve. With reference to the ASE 
curve (solid-blue curve in Fig. [2b]) obtained with ([3]) as a LOS probability function, the the value 
of the parameter a turns to be 1.15 within the range of A 1-50 BSs/km 2 , 0.48 within the range 
50-500 BSs/km 2 and 0.81 within the range 500-10000 BSs/km 2 . 


B. Frequency reuse 

To have a comprehensive view of the frequency reuse as an interference mitigation scheme, 
we need to assess the trade-off between the ASE and the network coverage probability, defined 
as 1 — P [7 < yth]. The results of this trade-off are shown in Fig. (3} where we plotted the network 
coverage against the ASE for different frequency reuse factors and base station densities. 

Firstly, we focus on the LOS/NLOS propagation; we can notice from this plot that, if we fix 
the BS density, higher frequency reuse factors enhance the network coverage but, on the other 
hand, determine a drop of the ASE. This is in line with what one would expect from frequency 
reuse. Nonetheless, if we have no constraint in the choice of the BS density, the ASE vs coverage 
trade-off improves as the frequency reuse factor increases. In fact, the trade-off curve we obtain 

4 If both serving BS and interfering BS are in NLOS with the user, the path-loss exponents of the serving BS-to-user channel 
and of the interfering BS-to-user channels are the same and, therefore, the power or the interference and of the received signal 
varies at the same slope as a function of the density. 
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(n 

2 


LU 
</> 
< 



—0— LOS/NLOS - Freq. reuse 1 
—0— LOS/NLOS - Freq. reuse 2 
—©— LOS/NLOS - Freq. reuse 3 
—©— LOS/NLOS - Freq. reuse 4 
-A- SL - Freq. reuse 1 
-A- SL - Freq. reuse 2 
-A- SL - Freq. reuse 3 
-A- SL - Freq. reuse 4 


0.65 0.7 0.75 0.8 0.85 0.9 

Coverage 


Fig. 3. ASE vs coverage trade-off for frequency reuse. The trade-off curves have been plotted for BS density equal to 10, 20, 
50, 100, 200, 500, 1000, 2000, 5000, 10000 BSs/km 2 , and compare the combined LOS/NLOS model with the single slope one. 

for a given reuse factor K lies on the top-right hand side with respect to the curve for reuse 
factor K — 1. This means that, by increasing the reuse factor and the base station density at 
the same time, it is possible to achieve better performance than with a lower frequency reuse 
factors; note, though, that this is true when there is no constraint in terms of BS density. 

By looking at the single slope PL curve in Fig. [3] it appears that higher frequency reuse factors 
should still be preferred in order to improve the ASE vs coverage trade-off. However, unlike 
with the LOS/NLOS path loss, increasing the BS density enhances the ASE with no loss in terms 
of network coverage. Yet, modelling the signal propagation with the combined LOS/NLOS path 
loss yields different results than with the single-slope PL. 

C. Partially loaded networks 

In this subsection we show the results of the cell densification for partially loaded networks 
with LOS/NLOS propagation. Differently from the case of fully loaded networks, we recall 
that a fraction of the BSs may be inactive and, thus, the density of interfering BSs Ai does 
not necessary follow the trend of BS density A (see Section IIII-EI and Fig. Q]). In Fig. [4a] and 
l4bl we show the outage probability and the ASE curves, respectively, as function of the base 
station density for difference user densities. To better understand the effect of the partial load 
on the network performance, we compare these curves with those for fully loaded networks. We 
reported the values of the probability p\ of a BS being active over the outage and ASE curves. 

We observe that, as long as pa > 0.9, the deviation from the fully loaded network case is 











Outage probability - P[y < y.] 
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(a) Outage vs base station density 



(b) Spectral efficiency vs base station density. 


Fig. 4. The probability pa given by 03 is reported as a dot on each curve. The outage probability has been obtained for 
7th = -8dB. 


minimal. However, as soon as A approaches the value of user density Au, the probability p\ drops 
and, as a consequence, the density of interfering Ai BSs grows slowly with A, up to the point 
where it saturates and converges to Au (see Fig. [Hi- At the same time, as A increases, the distance 
from UE to the serving BS tends to decrease, leading to an increment of the received power. 
Overall, the fact that Ai saturates whereas the received power keeps growing as A increases has 
a positive impact on the SIR; as a result, the outage probability (see Fig. [4a]) and the spectral 
efficiency improve once the density A approaches or overcomes Au- 

In regards to the ASE trend, we show the results in Fig. [4b[ According to (l20l) . the ASE trend 
is the combined outcome of the increase of the spectral efficiency and of the density of the 
active base stations. As the density of base stations increases and approaches the user density 
Au, the density of active base stations will converge to Au (see Fig. [4a]); given that the density 
of active BSs remains constant, the only contribution to the ASE increase will be given by the 
spectral efficiency improvement. As a matter of fact, we can see that, with respect to the case 
of fully loaded networks, the ASE curves show a lower gain when the density A approaches Au- 
To assess steepness of the ASE, we applied linear regression to the ASE curves in order to 
obtain the value of the parameter a corresponding to different intervals of A; we specifically 
consider the approximation for the curve corresponding to Au = lOOOUEs/km 2 (red curve in 
Fig. I4b1) . These values are a = 1.15 within the density range 1-50 BSs/km 2 , a = 0.43 within 
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Fig. 5. Transmit power per BS. The power has been obtained with an SINR threshold 7 th = —8dB and for tolerance Ag = 0.1%. 
The plot compares the TX power per BS for single slope-slope and LOS/NLOS path-loss for fully loaded networks. It also 
provides the curve for frequency reuse factor 2 and 3 and for partially loaded network with Au = 1000UEs/km 2 . 

the density range 50-500 BSs/km 2 and a = 0.46 within the density range 500-10000 BSs/km 2 . 
D. Transmit power per base station 

In Fig. |5] we show the simulation results of the transmit power per base station Ptx(A), which 
has been computed by using Algorithm [Has explained in Section ITV-AI In this figure we compare 
the results we obtained using the single slope and the combined LOS/NLOS path loss models. 

As we can see from this plot, the behaviour of the transmit power as a function of the BS 
density A is different in the two cases of single slope and combined LOS/NLOS propagation. 
With reference to Fig. 0 with single slope path loss, the power decreases linearly (in logarithmic 
scale) with the density; in the case of combined LOS/NLOS propagation, the transmit power 
exhibits different slopes as the as the base station density increases. We used linear regression 
to assess the slopes of the TX power curves (indicated by 5, as explained in Section IIV-CI) 
within different density intervals. With reference to the curve corresponding to fully loaded 
networks with LOS/NLOS propagation (solid-blue curve in Fig. [5]), the values of (Ft, 5) are 
(9.3-10 -9 , —1.9) within the A range 1-60 BSs/km 2 , (4.4 • 10~ 17 , —3.9) within the A range 60-300 
BSs/km 2 and (1.15 • 10 -9 , —1.44) within the range 300- 10000BSs/km 2 . 

The fact that the transmit power per base station decays more or less steeply with the density 
A depends on how quickly the interference power increases or decreases with A. As we explained 
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in Section ITV-Al the transmit power per base station -Ptx(A) has to be set so that the network is 
interference limited. Thus, if the channel attenuation between the interferer and the user decreases 
quickly as the density increases, a lower transmit power will be enough to guarantee that the 
interference power is greater than the noise power. In other words, if the interferer-to-user channel 
attenuation tends to decrease quickly as the density increases, so does the transmit power and 
vice-versa. For instance, for A G [60, 300]BSs/km 2 , the probability of having interferes in LOS 
with the user rises and, as a consequence, we have a lower attenuation of the channel between 
the interfering base station and the user. Hence, the / ; - rx (A) which guarantees the interference- 
limited regime will also decrease steeply with 5 = —3.9 as A increases. On the contrary, for 
A > 300 BSs/km 2 , most of the interferers will have already entered the LOS zone, meaning that 
the interferer-to-user channel attenuation drops less rapidly than for A < 300 BSs/km 2 ; for this 
reason, also -Ptx(A) will decrease less rapidly with <5 = —1.44. 

Let us note that, with increasing reuse factors N, the TX power decreases, as indeed a smaller 
bandwidth is used and, thus, the noise power is lower. 

E. Energy efficiency 

One of the most surprising outcomes of our study on LOS/NLOS propagation for ultra-dense 
networks is the effect of cell-densification on the energy efficiency of the fully loaded network, 
of which we show the results in Fig. [6a] The difference between the energy efficiency with 
single-slope and with LOS/NLOS path-loss is noticeable. In the case of single-slope PL, due to 
the linear growth of the ASE, t/ee (A) is a monotonically increasing function of the density A 
(see Section IIV-C1I) . In the case of LOS/NLOS propagation, from Fig. [6a] we observe that the 
energy efficiency exhibits a maximum, which is achieved for a given density A 0 . 

To explain this, we consider the case frequency reuse N = 1 (solid-blue curve in Fig. I6al) : 
from ([29b and with the values of the parameters P 0 (given in Table [[]), P T and i) (given in Section 
IV-DI) . and a (given in Section IV-AI) . the optimal point Ao is approximately 100BSs/km 2 . Beyond 
this point, the ASE gain is too low to compensate power consumption increase in the network, 
leading to a drop in terms of energy efficiency. From Fig. [6a] we can note that frequency reuse 
reduces the energy efficiency compared to N = 1. As a result of the lower ASE achieved at 
higher frequency reuse factors N, the energy efficiency drops as N increases. 

In Fig. [6b] we show the energy efficiency for partially loaded networks, for a user density Au 
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BS density [n.BS/km 2 ] 


(a) Energy efficiency for fully loaded networks. 



Fig. 6. (a) Energy efficiency vs BS density for fully loaded networks. The plot compares the energy efficiency for LOS/NLOS 

with single slope path-loss. The energy efficiency is given also for frequency reuse factors 2 and 3. (b) Energy efficiency vs BS 
density for partially loaded networks. Curves are given for Au = 1000UEs/km 2 and for three values of p. 

of 1000 UEs/km 2 . As we are dealing with partially loaded networks, we are interested in the 
BS densities A > Au, where energy efficiency strongly depends on the power saving factor p of 
the BSs in stand-by state. This is because the parameter p determines the energy saving of the 
inactive BSs, which become more numerous as the density A increases. Depending on the value 
of p, according to (1331) a local maximum may even occur at A* = Q p( U 1 _~^ • 

With p = 0.1 and with the values of a given in Section IV-Cl the local maximum turns to 
be A* = 7300BSs/km 2 . For higher values of p, X* is smaller than or too close to Au to be 
considered as a reliable estimate of a maximum; we recall from Section IIV-DI that this estimate 
can be reckoned as reliable only if A* is sufficiently greater than Au- In fact, we observe from 
Fig. [60 that there is no local maximum beyond Au for p = 0.3 or 0.6. 


VI. Conclusions 

In this paper, we have proposed a stochastic geometry-based framework to model the outage 
probability, the Area Spectral Efficiency (ASE) of fully loaded and partially loaded Ultra-Dense 
Networks (UDNs), where the signal propagation accounts for LOS and NLOS components. We 
also studied the energy efficiency of UDNs resulting from this propagation model. 

As the main findings of our work, we have shown that, with LOS/NLOS propagation, massive 
cell densification determines a deterioration of the network coverage at high cell densities, if the 
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network is fully loaded. Moreover, the ASE grows less steeply than a linear function at high 
cell densities, which implies that a larger number of base stations would be required to achieve 
a given throughput target with respect to the case of single slope path-loss. In regards to the 
energy efficiency, cell densification turns out to be inefficient for the network from an energetic 
point of view. In partially loaded networks, when the base station density exceeds that of the 
users, cell densification results in a coverage improvement. Overall, based on our findings, we 
can conclude that UDNs are likely face coverage issues in highly crowded environments with 
many users, which represent the worst case scenario for ultra-dense networks. 

Appendix A 

PDF OF THE DISTANCE TO THE SERVING BS 


Once the LOS probability function is known, from (fTTh we obtain the PDF of the distance to 
the closest BS as follows: 

P [r > R] = exp f - A f p L (\\x\\)dx\exp f - X f (1 - p L (||x||)) dxV (A.l) 

V Jb(P,R) J \ JB(o,d^(R)) J 

Assuming the integrals in (IA.1I) can be solved in a closed-form, with some symbolic manipulation, 
(IA.1I) solves in its general form as follows: 
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By taking the derivative of (IA.2I) . we obtain: 
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The PDF of the distance to the serving BS can finally be obtained as 
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UR) = [P [r > R]] = -P [r > R] Y f'M- 


(A.3) 


(A.4) 


m= 1 


If we assume the LOS probability to be given by ©, we can further develop (IA.1I) by solving 
the integrals in (IA.1I) and, with further symbolic manipulation, we obtain: 


_R? _ R eq 

P [r > R] = e nXL2e ^ ■ e ~ 7rXL2e ^ 
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where R cq = cZ eq ( R). Let us define the functions fi(R), HR), h (R) and their first derivatives 
f[(R ), f' 2 (R), and f 3 (R), respectively, as follows: 


fi(R) = nXL 2 e~^, / 2 (R) = -nX L 2 e4, / 3 (R) = -nXR 2 cq , f[(R) = -2ttA Re 


R* 
' L? 


f'(R) = vrA K 2 2f3 eq R 2 ^~ l e 


-K^R 2 Pe<i 


/'(fl) = -nXK&^R 2 ^-' 


By plugging (IA.5I) and f[(R), /^(-R), and / 3 (R) in (IA.4I) . we obtain the PDF of the distance to 
the serving BS. 

When the LOS probability function is given by ©, we obtain the PDF of distance to the 
closest BS station as follows. First, by solving the integrals in (IA.1I) and by some additional 
algebraic operations, we obtain P [r > R] as follows: 

P [r > R] = e 2n) ' L2e ~ T . e ^XLRe~T _ e ~n\R^ _ e -2nXL 2 e T? _ e ~2nXLR eq e E 3 - ^ g) 


Then, we define the functions /i(R), / 2 (R) • • • , /s(R) and we compute their respective deriva¬ 
tives f[(R), /a(R) ■ ■ ■ , /s(R) as follows: 

/i(R) = 2nXL 2 e~%, f[(R) = -2ir\Le~%, / 2 (R) = 2nXLRe~^, f 2 {R) = -2vrA(R-R)e~?, 

h{R) = —vrAR 2 q , /'(R) = -7rAR e 2 q 2/3 eq R 2Aq - 1 , / 4 (R) = -2 ttA 
/i(R) = 27rALR eq & q R^ q e-^^, / 5 (R) = -2 ttA LR cq e~^, 

K R^ e< i 

HR) = 2nXLK cq p eq R^- 1 (K cq R^ - L)e~ 

Finally, the PDF can be obtained by plugging f[(R), / 2 (R) • • ■ , /s(R) and (IA.6I) in (IA.4I) . 
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